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Abstract: Using gauge/gravity duality techniques, we discuss the sound-channel retarded 
rS correlators of vector and tensor conserved currents in a class of (2 + l)-diinensional strongly- 

coupled field theories at zero temperature and finite charge density, assumed to be holo- 
graphically dual to the extremal Reissner-Nordstrom AdS4 black hole. Using a combination 
fvq of analytical and numerical methods, we determine the quasinormal mode spectrum at fi- 

K^ nite momentum for the coupled gravitational and electromagnetic perturbations, and discuss 

f~^ the appropriate choice of gauge-invariant variables (master fields) in order for the black hole 

^^ quasinormal frequencies to reproduce the field theory spectrum. We discuss the role of the 

near horizon AdS2 geometry in determining the low-frequency behavior of retarded corre- 
f— V lators in the boundary theory, and comment on the emergence of criticality in the IR. In 

^^ addition, we establish the existence of a sound mode at zero temperature and compute the 

speed of sound and sound attenuation constant numerically, obtaining a result consistent with 
the expectations from the zero temperature limit of hydrodynamics. The dispersion relation 
rS of higher resonances is also investigated. 
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1. Introduction 

The recent surge in applications of the AdS/CFT correspondence in modehng strongly- 
interacting theories in the same universality class of systems relevant to condensed matter 
physics is an indication of gauge/gravity duality techniques being a promising new approach 
to study these systems.^ One particular such avenue, which we focus on in this paper, uses 
phenomenological classical gravity plus matter systems in order to study some aspects of the 
quantum critical regions in the phase diagram of condensed matter systems via holography. 
At a quantum critical point, the system undergoes a quantum phase transition at zero tem- 
perature and it is often modelled by a strongly-coupled conformal field theory, a regime where 
a standard field theoretic understanding is generically lacking. One may then hope that the 
AdS/CFT correspondence will be useful in addressing questions that go beyond the scope of 
the usual perturbative techniques. 

In this paper, we focus on a class of strongly-coupled (2 + l)-dimensional field theories 
at zero temperature and finite U(l) charge density. The simplest gravity background that 
captures the essential features of such theories in a holographic setup is the extremal Reissner- 
Nordstrom AdS4 (RN-AdS4) black hole. By studying the coupled electromagnetic and gravi- 
tational fluctuations about the RN-AdS4 black hole, the AdS/CFT correspondence allows us 
to calculate the retarded correlators of the U(l) charge current ( J^) and energy- momentum 
(T^iy) operators of the dual fleld theory for non- vanishing frequency and momentum. This 
holographic system has recently attracted much attention in the literature, inasmuch as the 
dual field theory exhibits a variety of emergent quantum critical phenomena, including non- 
Fermi liquid type behavior in the case of fermionic operators [5]. The emergence of the critical 
features in the IR physics of the dual theory is rooted in the fact that the near horizon region 
of the extremal background geometry is AdS2 x R^, and the assumption that there exists a 
so-called IR CFT dual to the AdS2 factor. 

As a complement to our previous studies [6, 7] of transport coefficients and transverse 
(shear) channel correlators of J'* and T^^ at finite frequency and momentum, here we ex- 
plore the structure of correlators of the vector and tensor currents in the longitudinal (sound) 
channel of the dual field theory. In particular, by using a combination of analytical and 
numerical techniques, we determine the spectrum of excitations for generic values of the (spa- 
tial) momentum. Along the way, we comment on the appropriate choice of gauge-invariant 
bulk variables (master fields), such that their associated quasinormal frequencies encode the 



'^See [1, 2, 3, 4] for accounts of applications in condensed matter physics. 



spectrum of excitations of the dual theory. Studying the retarded correlators at finite fre- 
quency and momentum, we present evidence for the existence of a branch cut on the negative 
imaginary axis, as well as a series of metastable modes corresponding to isolated poles in 
the lower-half complex frequency plane. At low frequencies (as compared to the chemical 
potential // conjugate to the charge density), but generic momentum, we find that the re- 
tarded correlators of the sound-channel vector and tensor operators exhibit emergent scaling 
behavior. As in the case of the shear-channel modes studied in [6, 7], these operators give 
rise in the IR to two sets of scalar operators with different conformal dimensions. 

Turning to the low frequency and momentum regime, we establish the existence of col- 
lective sound modes at zero temperature, and compute the associated speed of sound and 
sound attenuation constant. We compare our results with the expectations from the zero- 
temperature limit of well-known hydrodynamical relations, and discuss the extent to which 
both results agree. Similarly, we study the dispersion relations of higher resonances, which 
are fundamentally different from that of the sound mode: their frequencies uj remain finite 
as the momentum k approaches zero, so that the spectrum for to <^ fx and k <^ fi is effec- 
tively dominated by the sound mode. We provide further evidence to support this picture 
by computing the sound poles residue and comparing the Green's functions obtained numer- 
ically with a truncated spectral representation including the sound modes only, observing an 
excellent agreement for values of the momenta in the range < k < fi. 

This paper is organized as follows. In section 2 we review some notions concerning 
linear response theory (LRT), hydrodynamics, and the sound modes, and comment on the 
applicability of these concepts to the zero-temperature scenario which is the main focus of 
this work. In section 3 we review the RN-AdS4 black hole background, its near horizon 
AdS2 X R^ geometry, and the basic properties of the dual field theory. We then present the 
linearized Einstein-Maxwell system and introduce the gauge-invariant combinations known 
as "master fields", which allow us to reduce the system to only two degrees of freedom 
satisfying decoupled second order differential equations. In section 4 we exhibit the asymptotic 
solution for the various bulk modes, and redefine the master fields in order to obtain a new 
set of gauge-invariant variables which is better suited to the calculation of the holographic 
Green's functions. We then present formal expressions for the complete matrix of retarded 
correlators in the longitudinal channel, and elucidate their low-energy behavior by using 
the method of matching of asymptotic expansions, finding emergent scaling behavior of the 
spectral functions for generic momentum. Using the decoupled equations introduced in section 
4, in section 5 we numerically compute the quasinormal modes of the extremal RN-AdS4 black 
hole, which by virtue of our choice of master fields coincide with the dual field theory spectrum 
at zero temperature. Turning to the small frequency and momentum regime, we establish 
the existence of sound modes at zero temperature, and compute the speed of sound and 
sound attenuation constant numerically. Similarly, we investigate the dispersion relation of 
higher resonances (overtones) and discuss the extent to which the sound poles dominate the 
spectrum in the IR by computing their residues. We conclude in section 6. 



2. Field theory generalities 

Even though we will study the field theory spectrum at generic, not necessarily sniall^, values 
of the frequency and momenta, in order to interpret our results it is instructive to review some 
notions concerning the hydrodynamic (sound) modes, and comment on the applicability of 
these concepts to the zero temperature case which concerns us here. 

2.1 Finite temperature: LRT and the sound modes 

Here we will discuss only a few basic facts that will be relevant to our exposition. Detailed 
discussions of hydrodynamics in the context of the AdS/CFT correspondence can be found 
in [8, 9], for example. From the field theory point of view, the hydrodynamical regime is 
an effective description valid at wavelengths which are much larger than the mean free path 
of individual particles in the system. Among the various excitations in the spectrum, a 
prominent role is played by the so-called hydrodynamic modes, whose dispersion relation 
satisfies i^hydro(^ — t- 0) — ?• 0. The importance of these modes becomes apparent in linear 
response theory (LRT), where the response of a field (/>(t, r) (assumed to be a scalar field for 
simplicity) to a perturbation by an external source j{t,r) is given by 

m, ^) = - j dt'dV Gnit - t', f- f')Jit', r') , (2.1) 

where G/j(i — t\f— r ) is the retarded two-point function given by 

GR{t - t',f- r') = -ie{t - t') ([</.(t,f),,/.(t',f')]) . (2.2) 

We will denote by {w*} the set of poles of Gi?(t — t' ,r — f ) in the complex frequency plane, 
which define the spectrum of excitations of the theory. For simplicity, we assume that these 
poles are simple, and that there are no additional non-analyticities. Going to momentum 
space, we can then represent the retarded two-point function as 

Gr{uj, k) = 2, '^^ — I" terms analytic m [u, k) , (2.3) 

_r 1 w — ujjik) 



uJiGi 



and the response of the field to the perturbation represented by j{t' ,f ) takes the form 

m,f))=ie{t)l^e'^'-'^ Y, eI"^(-^)*-^I^«(-^)*Res(a;„fe)j(a;,,fc). (2.4) 



6{^*} 



Thus, l!ie{u)j) and lni{ujj) determine the frequency and damping, respectively, of the mode as- 
sociated with the eigenfrequency dj. In particular, eigenfrequencies with Im(;^j) > produce 



As compared to the chemical potential /x. 



instabilities that signal the breakdown of the linear response approximation. The eigenmodes 
with Ini(cjj) < represent the dissipative response of the system to the external perturbation 
and they are referred to as quasinormal modes (QNMs). Similarly, the residue Res(ti;j, k) de- 
termines the weight of the corresponding contribution to the total response. The important 
observation is that for small values of the momenta the hydrodynamic modes have, by defini- 
tion, the smallest imaginary part and therefore they dominate the relaxation of the system. In 
fact, for sufficiently large time scales^ the total response can be accurately approximated by 
the contribution of the hydrodynamic modes only. In particular, in the longitudinal (sound) 
channel one finds two hydrodynamic modes. Assuming rotational invariance in the spatial 
r = {x,y) plane, we can align the momentum with the rc-axis, i.e. k = kx; if the underlying 
theory has unbroken parity (y — t- —y) invariance, these two poles have the same imaginary 
part and opposite sign real parts, and obey a dispersion relation of the form 

Usik) = Csk-irsk'^ + 0{k^), (2.5) 

at small momentum, where Cg is the speed of sound and Tg is the sound attenuation constant. 
Denoting the residue associated with cOg by Rs(A;), the representation (2.3) for the correlator 
in the hydrodynamic regime can then be approximated by 

G1?""''(a;, A; ^ 0) ~ ^^^tt + ^^ + terms analytic in (w, k) . (2.6) 

u — ujs[k) uj + uj*[k) 

As one increases the momentum the effects of higher resonances become important, and the 
expression above should be modified accordingly. 

2.2 Zero temperature 

While it is tempting to use the familiar language introduced above, one should keep in mind 
that the main focus in this paper is the zero temperature scenario, where it is not clear in 
what precise sense one can talk about hydrodynamics. In this light, our results should be 
primarily understood as an exploration of the structure of retarded correlators of conserved 
currents in the longitudinal channel of a (2 + l)-dimensional strongly-coupled field theory at 
zero temperature and finite charge density. As discussed in [5, 6, 7], additional care has to 
be exercised when studying the zero temperature dual of the extremal RN-AdS4 black hole. 
This is due to the fact that the T — ?• and a; — )• limits do not, in general, commute. This 
implies that the naive T — t- limit of some finite temperature expressions does not agree 
with the corresponding quantities computed by setting T = initially. Additionally, the 
analytic structure of retarded correlators becomes richer; for example, isolated poles at zero 
temperature can now coalesce into branch cuts. 

While being cognizant of the subtleties mentioned above, in later sections we will find 
that some of our results do in fact agree with the T — )■ limit of certain hydrodynamic 



^See [9] for a discussion of the relevant time scale and the importance of the quasinormal mode residues in 
determining it. 



equations. In particular, at T = the notion of "hydrodynaniic limit" entails considering 
frequencies and momenta which are small compared to the scale set by the chemical potential 
/x. In this regime, we establish the existence of exactly two modes with a dispersion relation 
of the form (2.5), with the constants Cg and Ts matching, within the numeric precision, those 
predicted by taking the T — >• limit of well-known hydrodynamic expressions. This justifies 
referring to these excitations in the zero-temperature theory as "sound modes" . Just like in 
the finite temperature case, we will see that these modes effectively dominate the spectrum 
at small frequency and momenta: the retarded Green's functions computed numerically via 
holographic techniques will be shown to be in excellent agreement with the approximation 
(2.6). In fact, the fitting is robust not only for k <^ fi, but also for values of the momenta of 
the order of the chemical potential. For larger values of the momenta it becomes necessary 
to include the effect of higher resonances in order for the truncated spectral representation of 
the correlator to be accurate. 

3. The gravity background and the hnearized theory 

Consider the Einstein-Maxwell action in 3 -|- 1 spacetime dimensions with a negative cosmo- 
logical constant A = — 3/L^ 

S = ir2 I d^^V^ {R-2A- L^F.^.F'^'') , (3.1) 

2k^ J 

where L is the curvature radius of AdS4. The background we are interested in is the Reissner- 
Nordstrom AdS4 (RN-AdS4) black hole 

ds^ = g^.dx'^dx'' = ^ [-f{r)dt^ + dx" + dy^] + -^^dr^ , (3.2) 

A = Jl-—\dt, (3.3) 

which is a solution of the equations of motion obtained from (3.1) with 

/M = 1-m(^)%0=(^)\ ,.5™. ,3,4) 

Here tq is the position of the horizon, given by the largest real root of /(rg) = 0, so that 
M = 1 -|- Q^ in our conventions. The black hole temperature takes the form 

(3-Q2), (3.5) 



47rL2 
while its entropy, charge and energy densities are given by 



respectively [10]. This background is holographically dual to a (2 + l)-dimensional strongly- 
coupled field theory at finite temperature T and finite charge density p. The entropy and 
energy densities of the dual field theory are given by s and e in (3.6), respectively. Similarly, 
the asymptotic value of the bulk gauge field At{r — )• oo) = fi is interpreted in the dual theory 
as the chemical potential for the (electric) charge density. Little is known about the details 
of the dual theory from the field theory perspective. On the other hand, using holography 
a lot has been learned (especially thermodynamical properties) regarding its strong-coupling 
behavior; see [10] and its citations. 

The background becomes extremal when Q^ = 3; for this value of the charge the tem- 
perature is zero, but the entropy density remains finite. Since the solution is invariant under 
changing the sign of At we can choose /i to be positive, so that in our conventions Q = \/3 
at extremality. In this paper we will mainly work in the extremal limit, and refer to the 
corresponding dual theory as the "boundary field theory". As we will discuss in section 4, 
interesting properties of the boundary field theory in the IR limit stem from the features of 
the extremal near horizon geometry, which we briefly review below. 

3.1 Near horizon geometry at extremality 

Although the black hole temperature vanishes at extremality, its horizon area remains finite, 
a fact whose dual interpretation is that the boundary theory has a finite entropy density at 
zero temperature. In the extremal limit, /(r) in the background metric (3.2) takes the form 

/(r) = 1-4^1% 3 P)', (3.7) 



; 



which has a double zero at the horizon, and can be approximated near that region (to the 
leading order in r — rg) by 

m^^ir-rof. (3.8) 

The near horizon geometry at extremality is AdS2xR^. To see the emergence of this geometry, 
first change the radial coordinate r to rj defined by 

r - ro = — . (3.9) 



6r/ 



There is then a scaling limit [5] in which 



^ { j^2 . j„2\ , '''O /j^^2 , j,,2\ A Q 



ds' = -^\^-dt^ + dr]^j+^(dx^ + dy^j, A=^dt. (3.10) 

The curvature radius of the AdS2 factor is L2 = L/vG- The radial coordinate is interpreted 
holographically as the renormalization scale of the dual field theory, and the near horizon 
region corresponds to its IR limit. This implies that the AdS2 x R^ geometry encodes the IR 
physics (w — )■ 0) of the boundary theory. 



3.2 Linearized equations of motion: master fields 

The holographic calculation of the vector current and energy-momentum tensor operators of 
the boundary field theory entails solving the linearized Einstein-Maxwell equations for the 
corresponding fluctuations. We first define 

gfiu = (jfiu + hf_„y , A^ = A^ + a^, (3-11) 

where g^^ and A^ denote the background metric and gauge field, respectively, and /i^j^ and 
a^ represent the fluctuations. We work in the so-called radial gauge, where 

ar = 0, hr^ = 0, (3.12) 

with I' = {t,x,y,r}. We proceed by Fourier transforming the fluctuations 

V(i, X, r) ~ e-*"*e*'^-' V(0 , «/.(*, x, r) ~ e-'^'e'''^a^{r) , (3.13) 

where, without loss of generality, we have used the rotation invariance in the (x, y) plane to 
set ky = {) and defined kx = k. The fluctuations split into decoupled sectors depending on 
whether they are even or odd with respect to parity, y — t- —y. Accordingly, hty, h^y, ay have 
odd parity while hu, htx, hxx, hyy, at, ax all have even parity. The analysis of the odd parity 
(shear) modes was performed in [7]. In this paper, we study the even parity modes which 
translate into the longitudinal modes of the boundary theory. 

It is convenient to raise indices in the metric fluctuations by the background metric ^^j/, 
and work instead with /i*t, /i^t, h^x, h^y. We also define the dimensionless quantities 

(3.14) 





r 




CO 




k 


u = 


1 


vo = 


~ 5 


q = 


= — 




ro 




^^ 




/^ 



The linearized Maxwell equations read' 
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= 2u^f{u) [u^a'^{u)\^ - 2Q2q [qai(ii) + xoax{u)] 



-fiu^fiu) [h'[{u) - h-'x{u) - h^'y{u)\ , (3.15) 

= u^f{u) [u^f{u)a'^{u)\' + Q'^vo [qat(n) + t»a^.(n)] + i^u^ f{u)h''[{u) , (3.16) 

= 2u^ [voa[{u) + q/(uX(n)] - fim [h\{u) - h'^xiu) - /i\(n)] + 2/iCi/i^t(u) , (3.17) 

where f{u) is given in (3.7). Notice that not all of the above equations are independent; 
equation (3.16) can be obtained from equations (3.15) and (3.17), for example. For the 



*We keep the black hole charge Q arbitrary (0 < Q < y/S) in writing the equations in this section. In later 
sections, when we discuss the extremal case, we set Q = v^- 



linearized Einstein equations we obtain 

= f{u)l2uy{u)h'';{u) + [lOnV(u) + 3nV'(n)] h'[iu) 

+ [u^fiu] + 2u'fiu)] [h^Uu) + hy'yiu)] I - 8fi-'Q\^f{u)a[{u) (3.18) 

+ 2g2/(w) (2 - q^n^) /i*i(u) + 2Q2tx,2n2 [/,a=^(^) + ^j/^(^)] + 4Q2qtx>u2/i^(n) , 



= fiu)lu%^:iu) + [bu'fiu) + u'f'iu)] r;(n) + u'fiu) h'',{u) + hy'y{u) 

Wu'h\{u)-Q^f{u){2- 
+ 2Q2qt»n2/i^i(u) + 4;U~^Q\2^(n)a;(n) , 



+ Q2«,2^2^^^.(^) _ g2^(^) (2 + q2^2^ ^t^^^^ _ Q2q2^2^(^)^j/^(^) (3_^9) 



= f{u)lu^hy';^{u) + [5nV(n) + nV'(n)] /^^'.(n) + u^f{u) [h'[{u) + /i-^(n) 

- 2Q^f{u)h\{u) + Q^u^ (nj2 _ /(u)q2) /^^/^(n) + 4/x-iQ\V(n)a;(n) , (3.20) 

= f{u) [u^h'^iiu) + Au^h%{u)\ + Q2qn,n2/i%(n) + A^-^Q'^u^ f{u)a'^{u) , (3.21) 

= 2u^f{u) \h't{u) + r'^(n) + hy'^{u)] + [uV(^i)]' ki(^) + h^'xiu) + hy'y{u) 
+ 2u^f'{u)h^[{u) + 4Q2/it^(u) - 8/i~^g2^2Q;(u) , (3.22) 

= 2vofiu) [/i";(u) + hy'yiu)] + 2<if{u)h%{u) - wf'iu) [h\{u) + /i\(7x)] (3.23) 

-2^f{u)hMu), 

= 2qu^f{u) \h^[{u) + hy'yiu)] - 2wu^h%iu) + qu'^f'{u)h\{u) 

- Sfi-^Q"^ [qat(n) + voa^{u)] . (3.24) 

The linearized Einstein-Maxwell equations presented above can be reduced to a set of 
two decoupled second-order differential equations for the so-called "master fields" ^±{u), 
introduced by Ishibashi and Kodama in [11]. In appendix A we explain in detail how to 
obtain the decoupled equations for the master fields. Basically, one defines the gauge- invariant 
quantities 



1 



A{u) = u a[{u) + fj,hyy{u) fih i(n), 

$(«) = uhyyiu) - [QW + u'fin)]-' nV(n) [h^M + /i%(tx)] , 



(3.25) 
(3.26) 



in terms of which the equations (3.15)-(3.24) reduce to a set of two coupled second-order 
differential equations given by (A. 3) and (A. 4). The resulting two equations can be decoupled 
by introducing the master fields^ 

^±iu) = a±{u)<^{u) + fi-^Q^Aiu) , (3.27) 

where we have defined 

a±{u) = \f±{Q) - ^ , F±(Q) = ^ [(1 + g2) ± V(l + Q')' + (16/9)Q^q2l • (3.28) 

2 U 4:1 J 

The decoupled equations for the master fields then take the form 

u^fiu) [tx2/(n)$'±(u)]' + [Q2n,2 _ U±{u)] $±(n) = 0, (3.29) 

with 

C/±(u) = ±7/(^) |a±(n)y(n) + -^ [QWH{u) + 8Q2/(^)] + — [Hiu) - 2Q^^] 
I 4:a±[u) ■' u '- -' 

+ ^ [48u'f{u)ifin) - 1) + 20txV(n)/'(n)] | , (3.30) 

where 7, V{u) and H{u) in (3.30) have been defined in appendix A. In the special case where 
the (electric) charge of the black hole background vanishes, namely Q = 0, the electromagnetic 
and gravitational perturbations decouple, and the master fields introduced above reduce to 
a single field carrying the information from the Einstein sector. In our notation, this field 
is given ^+{u) = $(m), and it satisfies equation (A. 19). The gauge sector is dealt with by 
introducing a second independent gauge-invariant mode. 

In the next section, we will see that the Ishibashi-Kodama master fields introduced in 
(3.27) are not well suited to the analysis of the sound-channel retarded correlators of the 
boundary theory via holography.^ As we explain below, this is because of the relationship 
between the leading term in the asymptotic expansion of ^±{u) and those of the metric and 
gauge field fluctuations. The guiding principle is that, in gauge/gravity duality, the ambiguity 
in the choice of asymptotic boundary conditions for the gauge-invariant fields is removed by 
demanding the quasinormal mode spectrum in the gravity side to coincide with the poles of 
the retarded correlators in the dual field theory [12]. 

The observation that the Ishibashi-Kodama-like master fields with asymptotic Dirichlet 
boundary conditions may not be appropriate for the analysis of the sound-channel correlators 
via holography was first made in [13], for the case of an uncharged AdS4 black hole. In order 
to obtain the sound mode in the boundary theory, the authors of [13] proposed that one 



^The master fields defined in (3.27) differ from the ones given in [11] by a multiplicative constant. 
^They are a convenient choice, however, to analyze the shear-channel correlators [7]. 
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should instead impose a Robin ( "mixed" ) boundary condition on the master field (recall that 
in the case where the background is uncharged, ^± reduce to a single field which carries the 
information about the gravity fluctuations). Our approach here is different: we introduce 
new master fields by taking linear combinations of $±(m) and their derivatives, such that the 
leading term in the asymptotic expansion of the new master fields contains the leading term 
in the asymptotic expansions of the metric and gauge field fiuctuations only (which are the 
sources for the corresponding operators in the boundary theory). In other words, we perform 
a canonical transformation in the space of master fields and their conjugate momenta, so that 
the new master fields satisfy a Dirichlet boundary condition asymptotically. The advantage 
of this approach is twofold: on one hand it makes the relation between the asymptotics of 
the master fields and the gauge- variant sources more transparent, while at the same time it 
allows us to use the various known techniques for the calculation of the quasinormal mode 
spectrum, such as the (Leaver's) matrix method we will employ later on. 

4. Green's functions and criticality 

In this section we present the asymptotic expansions of various bulk fields and the renormal- 
ized boundary action, and use them to obtain formal expressions for the retarded two-point 
functions of the vector current and energy-momentum tensor operators of the boundary the- 
ory in the sound channel. We then analyze these expressions in the small frequency limit and 
comment on the emergent critical behavior stemming from the properties of the near horizon 
geometry of the extremal background. 

4.1 Asymptotic expansion of the bulk modes 

Asymptotically as u — )• oo, the solution of the Einstein-Maxwell equations (3.15)-(3.24) have 
the expansions 

a^(n ^ oo) = a^ [1 + O (u^^)] + ^vr^ [l + O (u'^)] , (4.1) 

hf'^iu ^ oo) = /i^^ [l + O {u^^)] + ^TT^"^ [l + O (n-2)] , (4.2) 

where it is understood that a^, vr^, /I'V and vr'^j^, hence a^{u) and h^u{u), are all functions 
of vo and q. We have denoted the leading constant terms in the asymptotic expansions of 
the modes by a hat. These are the sources for the corresponding operators in the boundary 
theory. Inserting (4.1) and (4.2) into the six second-order Einstein-Maxwell equations, we 
find 

1 g2 q2 



1 

^24^ 



Gvr^ + 6(1 + Q2) t,^ + q4jj, ^^2 _ q2^ (q^^ ^ ^^^^ ]+•••, (4.3) 
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at{u — )• c«) 



1 



Q' 



Q' 



°* + n^*+2^^('^* + "'"-^ + 24^3 



4q(q7ri + u>7ra; - qh\) 



.2?s 



+ 



1 



q2/i*^ + (a>2 



c^) n 



24^4 



3 K"x + ^\ 



+ QS (q^ - O)^) (qa* + WQi 



+ 



(4.4) 



02 r ^ 

hMu ^oo) = h\ + ^ 2tx)q/i^4 + tx)2/i^^ - q^h\ + (txj2 + q2) /j?/^ 



4m2 



+ ^vr'i 



161*4 



247ri + 12h\ - QW (2tx)q/i^j + tw^/i^^ - q'^h\ - {w^ - q^) /jJ/^] 1 + . . . , (4.5) 



/i^^(u ^ oo) = /i^^. + ^ r2njq/i^j + nj^/j^^ _ Cf^h\ - (rx)^ + q^) ^ ^ ^ 



4n2 _■ 






+ 



16^4 L' 



SvTi + 4/i*i - Q2«^2 ('2tx)qA^i + tx)2/i^^ - q2/^*^ _ (^2 _ ^2^ f^y\]^...^ (45) 



/i^y(M — ;> 00) = /i^, 



Q 



2 r 



?/ 4^2 



2tx)q/i^4 + n)"/i^^ - q2/i'j - n?" - q" /i?^ 



_i -jT-y 



Q' 



16«4 



2tx ^2lt 



87rt+Ah\ + Q' (w' - q") 2«)q/i^, + n)^^^ - <f'h\ - (xx>' - q" /i^, 



+ 



(4.7) 






+ 



16^4 



16^:,. + g2tx,q (2tx)q/i^i+tx)2/i^^-q2/i*^_ (t»2-q2)/iS/J +.... (4.8) 



When expanded asymptotically, the four constraint equations (the uu, ut, and wx-components 
of the Einstein equations together with the u-component of the Maxwell equations) yield 



4Q2 



qvr^ + qvr^t - twvr^ = TT (l + Q^) ^h\ 



{<\at + voax) 



q7r''t+t»7r'''^. + nj7r^. 



7r*i + Ti\ + ^% = . 
Using (3.29), the asymptotic expansion of the master fields ^± takes the form 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



$±(m^oo) = ^1°N 1 + 
+ ^|l + — 

u 1 6ti2 



1 

2^ 

8F2 
g4q4 






i^± 



3Q2q2^3 

'4 



16F| 



2^2 



6F±u3 



16Ff 



6„6 



Q6q 



Q4q4 



+ Q'q 



+ 



+ •• 



(4.13) 
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From the definition of tlie master fields in (3.27) and tlie asymptotic expansions of at{u), 
h*t{u), h^x{u) and hyy{u) given in (4.4), (4.5), (4.6) and (4.7), respectively, one concludes 
that 



$ 



(0) 



1 



QH 



2r,2 



1 



F±tt\ - li-'QWT^t + Fih\ - -QWh\ 



$ 



(1) 



■^ {-lF±-'t - f^-'QVn + Fihyy - lQVh\ 



1 



- 7^±^/i - fJ' Q z, 



(4.14) 



(4.15) 



where we have defined 

Za = Q^q (qa* +t»aj 



Zh = Q" \lX0<^h\ + tX>2/i-^ - q2/,t^ + (q2 _ ^2^ f^y^ 



(4.16) 
(4.17) 



Notice the appearance of the momenta 7r*t and lit in <1>^ : this is the reason why the Ishibashi- 
Kodama master fields ^±{u) are not appropriate for the holographic analysis of the boundary 
theory Green's functions in the sound channel. As we mentioned before, this issue is specific 
to the sound channel and does not arise in the shear channel. 



4.2 New master fields 

From equations (4.14) and (4.15) we obtain 



|^*l°> + *i" = -lF,4-M-'Q=5„. 



(418) 



which is independent of the momenta 7r*i and vr^. This suggests that we define two new master 
fields 



^±{u) = ^^±{n) - u^f{u)^'^{u) , 



Q2q 



(4.19) 



whose leading asymptotic (constant) terms do not depend on the field momenta. Using the 
fact that ^±{u) satisfy (3.29), the decoupled differential equations for ^±{u) take the form 



where 



[u^fiu)^'^iu)]' + J±iu)^'^{u) + K±{u)^±iu) = , 



J±iu) 
K±{u) 



4F| 



-I -1 






u'f{u)U'^iu) 



v?f{u) 






(4.20) 

(4.21) 
(4.22) 



13 



Asymptotically, ^±(m) take the form 



*±(n ^ oo) = ^I/f [l + O (n-2)] + ivI/W [l + O (n-i)] . 



(4.23) 



Using (4.19), (4.14) and (4.15), the two coefficients ^j_ and ^j_ are related to the asymptotic 
values (sources and momenta) of the metric and gauge field fluctuations through 



^>± ^ mf 



* 



± 



1 

'i 
3 

2F± (I 



(1) _ c± 



F±Zh - fJ- ^Q^Za , 



A^2lt 



F±7t\ - fi^'QS'TTt + Fihyy - -Q\'h 



Q2q2 \^4 



F±Zh + Ai Q Za 



(4.24) 



(4.25) 



where Za and Zh are given by (4.16) and (4.17), and we have defined 



Arp2 



(4.26) 



In order to compute retarded correlators, the Lorentzian AdS/CFT recipe of [14] instructs 
us to apply infalling boundary conditions for the bulk fields at the black hole horizon. Im- 
posing infalling boundary conditions for '^±{u) at the horizon, ^j_ will be related to '^±. 
Parameterizing this relationship as 



* 



(1) 

± 






^±, 



(4.27) 



equations (4.24) and (4.25) reduce to 



Q'q'n±^± = -\f^7^\ - f^-'QS^7Tt - lQVh\ + F^hyy 



(4.28) 



Solving the above two equation for 7r*j and ttj, we obtain 



vr*t = (1 + g2) hVy + 



1 Q'q 



2„2 



6 F+ - F„ 



(F+n+ - F_n_) 4 + ifi-'Q' (n+ - n_) za\ , (4.29) 
fiQW (n+ - n_) 4-4 (F_n+ - F+n_) d . (4.30) 



II-i- will be determined numerically. The other four field momenta are easily obtained by 
substituting equations (4.29) and (4.30) into equations (4.9)-(4.12). 
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4.3 Boundary action and retarded Green's functions 

The renormalized action in our case is given by [15] 

6 



Z/^4 Of, 



d'x^/^g [R+j2- L^P^^P""" 



[ Sx v1^2i^ - y / <fx^\-L I Sx^M ^^^R , (4.31 

JdM ^ JdM JdM 



where 7^1/, K and ^^'R are, respectively, the induced metric, the trace of the second fundamen- 
tal form and the intrinsic curvature of the 3-dimensional boundary dM 7 The counterterm 
proportional to ^^'R is introduced in order to remove a linear divergence in the on-shell ac- 
tion for the RN-AdS4 background. To obtain Green's functions holographically, it suffices 
to consider the (gravity) action up to quadratic terms in the fluctuations. The renormalized 
on-shell action then takes the form 



c _ 0(1) _|_ 0(2) 



(4.32) 



where the S^eL piece is linear in fluctuations and determines the one-point functions of the 
dual operators, while Srcn is quadratic in fluctuations and determines the two-point functions 
of the conserved currents in the dual field theory. In our case, ^ren is worked out to be 



sill = ' / (fx {T\)Uj + / (fx {J') at 



dM 



(4.33) 



dM 



with the one-point functions given by 



(n) 



lim — 



n 



3 



« + iM', 



2r2 

(J*)= lim ^\[duAt) 



(4.34) 
(4.35) 



where the bar denotes quantities that are evaluated in the background and i,j = t,x,y. In 

(2) 
the same way, Sren is given by 



S^'^l = C lim / d^x 
"^°° JdM 



u^f 



+ ^h''th% + 4^n2 {ata't - fa^a'^) + £cont \ , (4.36) 



''Tf ^A" 



If n^ denotes the components of the outward pointing unit normal vector to the (timeUke) boundary, we 



have 7^^ = g^^ - n^n„ and K^^, = -(l/2)£„7p^ = -7'',tVpn^ 



K = T^'J^M- = -V^n" 
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where C = rQ/{4:K\L ) and the "contact term Lagrangian" £cont is given by 



C 



u 



cont 



VJ'^ ' V ^^ 



dxdt^.^ 52 _ , ^52 ,\ O 



Q^ — /i% 2^h^t - -^h% + f^h 



/i- 



fi- 



^ 



a^h^'t + at ht 



n'f 



1 



1 iif'\ 
tr^i ' -I ^ ' ^ \h\ {h% + hyy) 



1 ^ 1 h\h\ + 2 -1 + ^ 



+ (1 - -^ + ^"j ih\h\ + hyyhyy) + 2 (-1 + 






4 / 


1 \ 1 




1+ ^ /l^/l^ 


/ 1 


^/7y 



(4.37) 



In momentum space we write 



^ 2n f ^ f dtodq 



5(2) 



tt\ h\ + T^yy h\ + tt\ h\ + T^yy h^^ + 7r*i hy. 



+ ^^^ hyy + 2^^i h'^t + -^ (vTi Oi - TT^ a^) + £contJ , (4.38) 



where >Ccont denotes the Fourier transform of the contact term Lagrangian, and it is under- 
stood that TT^x h^t = ^^x{'<^,<\) ^*t(— 1^) —'^) and so forth. Defining 



Gi(tx>,q) = C 
G2(i»,q) ■■ 



F+-F. 
8C Q^ 

~F:-F^ 



(tX)2-q2)2(F+n+-F_n_) + 3C7(l + Q^)(2 ^ 



w 



fJ- 



(t»^ - q^) (F„n+ - F+n_ 



(4.39) 
(4.40) 



we can write the various retarded Green's functions in a compact way. For notational simphc- 
ity, the contact term contributions are not displayed in the expressions below; their explicit 
form is given in appendix B. Using the Lorentzian AdS/CFT recipe of [14, 16], for the 
retarded two-point functions of the form {Tfj^uTa/s) we find 



1 q^ 
Gii,ti(t»,q) = - ^ ^_ ^ Gi(to,q) , Gxx,tt{^^'^) 



1 q2t»2 



Gyy,tt(t»,q) 



(-Jxx,xx{^i^) 



2(t»2-q2)5 
1 q' 



2tx)2-q2 



1 tx>^ 



2(tX>2-q2)2 



Gi(t»,q), Gx.t,tt(a>,q) 



Gi(rx>,q), Gyy^xx{^,<f) 



2(tX)2-q2)2 

1 q'^o) 



2(t»2_q2)2 

1 r»^ 



Gi(tx,,q), (4.41) 



Gi(«>,q), (4.42) 



2n?2_q2 



Gi(tx>,q), (4.43) 



1 qo?"^ 1 

G^t,xx(t»,q) = --T^y— ^Gi(t»,q), Gyy,yy{vo,c^) = -Gi{m,q), (AAA) 

G'xt,yj/(t»,q) = -^ :,Gi{m,c\), Gx-t,xt(t»,q) = ^^ ^-^Gi{xx>,q). (4.45) 



2nj2_q2 



(t»2 — q2)2 
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Similarly, the current-current correlators of the form (J^Ju) read 



Gt,t{r«,<f)= / 2 ^2(«^,q), (4.46) 

G'x,x(n',q) = ^ 2<^2(t»,q), (4.47) 

t»^ — q^ 

G,,,(t»,q) = --^^G2(tx),q). (4.48) 



Defining 



G^ix(j»,<l) 



4Q2 



2 ^2^ 



-Gi(a),q) ^ ^G2(ro,q, 

/x 8 



(4.49) 



3(1 + g2) 

(recall that we set Q^ = 3 at extremal! ty) the "mixed" correlators (J^T^p) are given by 

q^ q'^tw 

Gt,it(t»,q) = -^^ ^v^G',„i,(a),q), Ga.,it(i»,<l) = ^^ ^v^G^ix(n',C|) , (4.50) 

(t»^— q^)^ (tx)^— q^j^ 

q'^tx)^ qtx>'^ 

Gt,xx{'^,'^) = --r^ ^T^G,„i,(tx),q), Gx.,a.x(i»,<l) = r^ ^T^Gn,i,(nj,q), (4.51) 

(t»^ — q^)^ (t»^— q^j^ 

q2 qn) 

Gt,TO(t»,q) = ^2^-^G^ix(t»,q) , G^,yy(t»,q) = 2^-^G^ix(a>,q) , (4.52) 

q'^tx) q^tx)^ 

Gt,xt(a',q) = T^ ^T^G^ix(n',q) , Gx,xt{'x>,q) = --j-^ ^T^G„,i,(rx>,q) . (4.53) 

(uj^ — q^)^ (tx)^ — q^j^ 

Naturally, the matrix of Green's functions is symmetric, namely Gxx,tt{^,'^) = Gtt,xx{^,'^), 
Gt^ti (tx),q) = Gij_i(txj,q) and so forth. The form of the correlators above is consistent with 
the field theory Ward identities (see [2, 12], for example), and it agrees with previously found 
expressions in the simpler case where there is no mixing between the electromagnetic and 
gravitational perturbations [17]. It is worth noticing that the apparent pole at tx? = itq in 
the expressions for the correlators is spurious, and reflects only the choice of normalization 
of Gi and G2. 

4.4 Small frequency expansion 

In this subsection we obtain analytical results for the small frequency behavior of the retarded 
Green's functions, by expanding ^-i-(tx),q) and n-i-(tx),q) for small tx>. We write ^±{u) as a 
power series in tx>, 

^o±{u) = ^gi(«) + xv^'iiiu) + t»2^gi(u) + • • • , (4.54) 

with the subscript denoting the outer region of the background geometry (see the discussion 
below). At extremality, extra care is needed in taking the limit tx) — )• near the horizon. 
This is due to the fact that in the extremal case f{u) has a double pole at the horizon. A 
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convenient way of handling the Hniit was introduced in [5]. Basicahy, one reahzes that near 
the horizon the equations (3.29) organize themselves as functions of C = "^^j where, using 
(3.9) and (3.14), 

u=l + -^ . (4.55) 

\/l2C 

Since the coordinate C is the suitable radial coordinate in the AdS2 factor of the near horizon 
region, the t» expansion of ^± in that domain is given by 

*/±(C) = ^S(C) + «)*2(C) + «)'^g(C) + • • • , (4.56) 

with the subscript / denoting that the fields ^± are expanded as a power series in to in the 
near horizon geometry (inner region). After imposing infalling boundary conditions for ^/-|- 
at the horizon, we match the inner and outer expansions in the so-called "matching region" 
where the C ~^ and t»/C — >• limits are taken. Since equations (3.29) are linear, if we 
require that the solutions for the higher order terms in the expansions (4.54) and (4.56) do 
not include terms proportional to the zeroth-order solutions near the matching region, we 
just need to match "iffliO to ^o±(m) [5]. 

4.4.1 Inner region 

Near the horizon, we find that ^j^{C) satisfy 



■S"(0+(5^±^^^^i±5-ll<i( 



-* T(C) + ^L^±i;vi^:L _ 1 * ".,0 = , (4^57) 



which are identical to the equations of motion one obtains for two massive scalar fields in 
AdS2, with masses m± given by 

2 

mlLJ = 1 + |- T 7rTq2 . (4.58) 

Note that L2 in (4.58) is the curvature radius of AdS2- 

Assuming that there exists a CFT dual to the AdS2 region, one can think of ^]'ji(C) 
as being dual to scalar operators 0± in the IR CFT. According to the standard AdS/CFT 
dictionary, the conformal dimensions 6± = u± + 7^ oi the operators 0± are related to the 
masses m± of the scalar fields ^j^ by 

i^± = ^^1 + 4m|^2 = 1 ^5 + 2q2 ^ 4^/1+^ . (4.59) 

Next, we impose infalling boundary conditions for ^j_|:(C) at ^ — )■ 00. Having done so, near 
the matching region (where one takes the C ~^ and tx>/C — )• limits) ^jj:(C) take the form 

^fliu ^ 1) = \{u - l)-i+'^± + g±{xx>){u - ir'2-^^] , (4.60) 
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where we used (4.55) to change C back to u. Notice that a specific normahzation for ^]^_|i 
has been assumed in writing (4.60). Clearly, such a choice does not affect the calculation of 
the boundary theory Green's functions. The functions Q±(rx>) in (4.60) denote the retarded 
Green's functions of the IR CFT operators 0±, and are given by [5, 6] 

r(l-i/±) /ix)\2;^± 



g_t(t») = -2i/±e" 



■nTi/± - 



r(i + i/±) V2 



(4.61) 



4.4.2 Outer region and matching 



,(0) 



In the outer region, the equations for ^^^(li) are obtained by setting rw = in (3.29). The 
solutions for ^Q^(n) are given by a linear combination of (u — 1)^2+^^* and {u — 1)~2~^± 
near the matching region. Let us define 

r,^°\u) = {u- l)-"2+-± + . . . , ^(0)(n) = (n - l)-"^-^ + • • • , 



u 



(4.62) 



Then, matching ^q^(u) to (4.60) we obtain 



,(0) 



.(0) 



(0), 



*0±M= vTi^)+0±ito)e±'iu 



Following the discussions in [5, 7], to higher orders we can write 
'n±{u) = r]^ {u) +vor]^_^'{u) +w'^r]]_'{u) -\ , 

e±(^) = ef («) + tx'd'V) + «>2d'^(n) + • • • , 

where t/^ (u) and ^^T" {u) are obtained demanding that in the u 
distinct from ry^. (ti) and ^j_ (w), respectively. Thus, 

^o±{u) = [v±{u) + a±(t»)e±(u)] . 



An) 



Near u — )• oo, one can expand ?/j_ (u) and ^^ (u) as follows (here n > 0) 



(n). 



(it.) I \ 

r]j_[u^ ooj 



,{") 



1 + 



u 



4-) (^ ^ 00) = 4") (1 + • • • ) + 4^^) - ( 1 + 



u 



(4.63) 

(4.64) 
(4.65) 

1 limit, they are 
(4.66) 

(4.67) 
(4.68) 



where the coefficients aj!l ,b^ ,c^ and d^ are all functions of q. Thus, asymptotically we 
obtain 



*o±(^i -^ 00) ~ ^± f 1 + — 1 , 



where 






,(0) 



,(1) 



M 



.(1) 



6f + «,6W + O (a,2)l + g^(i„) [40) + u,^!) + (,„2y 



(4.69) 

(4.70) 
(4.71) 
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4.5 Criticality: emergent IR scaling 

In general, the coefficients a}^ ,b^]^, c^, d^ should be computed numerically. Once these 
coefficients are known, one can substitute the expressions for ^± and ^±11-1- into (4.41)- 
(4.53) in order to compute the retarded Green's functions. In the small frequency (but 
generic momentum) regime, the Green's functions can be analyzed semi-analytically even 
though the coefficients o^ , b^ , c^ , d^ are not known explicitly. Since the coefficients 
a)^ ,b]^ ,c^ , d^ are all real, the complex parts of ^-t and ^±II± are determined by the 
IR CFT Green's functions ^-|-(t»). Assuming that the product a_^_ a_ does not vanish for 
generic values of momentum q, expanding the denominator ^_|_^_ in (4.39) and (4.40), and 
noticing that 2i/_ > 3 while 2z^_|_ > 1, we obtain the following small frequency scaling behavior 
for Im Gi(tx),q) and Im G2(ix),q): 



Im Gi(tx>,q) = q2 ( 1 + ^/l+^) eo(q) Im g+{vo) [! + •••]« xx>^''+ , (4.72) 



Im G2(tx>,q) = ^ (l - \/l+q2J eo(q) Im g+{vo) [! + •••]« t»2^+, (4.73) 

where the dots represent subleading terms which vanish as tx> — ?• 0, and 

As a result, the matrix of the spectral functions of the boundary theory operators Jt, Jx, Tu, 
Txx, Tyy, Txt shows scaling behavior at small frequency 

Im Gtt,tt oc t»2^+ , Im Gxx,tt oc t»2+2!'+ ^ im Gyy^tt (xxo^''+, 

Im Ga;i,tf oc t»i+2'"+, lm.Gxx,xxOcvo'^+'^''-, Im Gj;y,^^ oc t»2+2'^+, 

Im Gxt,xx oc t»3+2i^+^ ini Gyy^yy OC t»2''+, Im Gxt,yy oc xo^+'^''+ , 

ImG,t,,iOct»2+2-+, ImGt,itoct»2-+, Im G,,it oc n>i+2-+, (4.75) 

Im Gt,xx oc nj2+2'^+ , i^ Gx,xx oc nj3+2'^+ ^ i^ Gt,yy oc tx>2'^+ , 

ImG,,yyOcnji+2.+ ^ Im Gi,,i oc nji+2-+, Im G,,,t oc t»2+2^'+, 

Im Gt,t oc tx)2'^+ , Im Gx,t oc xo^+'^''+ , Im Gx,x oc vo'^+^''+ . 

Notice that the low-frequency scaling behavior of the above spectral functions is emergent: 
the scaling is a consequence of the near horizon geometry (which translates into the IR physics 
of the boundary theory) containing an AdS2 factor. Although eo depends on details of the 
background geometry in UV, the scaling behavior rx>'^'^+ does not, as long as the AdS2 part 
of the geometry is unchanged. 

5. The spectrum of the boundary theory 

As we have discussed, the set of poles of the Green's functions (4.41)-(4.53) in the com- 
plex frequency plane, hereafter denoted by t»^(q), are generically given by the solution of 
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'I'-l-(t»,q) = (with infalling boundary conditions at the horizon). Thus, we have set up 
the problem in such a way that finding the spectrum of the boundary theory amounts to 
finding the quasinormal frequencies of ^± in the RN-AdS4 black hole background: in the 
context of the AdS/CFT correspondence, this should be regarded as the definition of the 
quasinormal frequencies in asymptotically AdS geometries [18, 19, 14, 12], thus removing the 
ambiguity associated with the choice of asymptotic boundary conditions which is inherent to 
the eigenvalue problem in the gravity side. 

Aside from its connection to the dual field theory retarded Green's functions, finding the 
electromagnetic and gravitational quasinormal frequencies of the extremal RN-AdS4 black 
hole is an interesting problem in its own right. From the gravity point of view, knowledge of 
the QNM spectra allows one to construct the dissipative response of the black hole geometry 
to small perturbations. In the case of the non-extremal RN-AdS4 black hole, such quasinormal 
frequencies were numerically computed in [20]. For the extremal RN-AdS4 black hole, the 
quasinormal frequencies of the odd-parity (electromagnetic and gravitational) modes (which 
are dual to the operators in the shear channel of the boundary theory) were recently computed 
in [7] . In this section, we extend our previous computation to the case where the perturbations 
are even under the y — )■ —y parity reflection. Using gauge/gravity duality, these fluctuations 
are mapped to operators in the sound channel of the dual field theory. See [21] for an 
extensive review of the subject of quasinormal modes in general relativity and the AdS/CFT 
correspondence. 

5.1 QNMs: matrix method 

Although there exist various semi-analytic methods for calculating small or large (magnitudes 
of) quasinormal frequencies, it is often necessary to perform numeric calculations in order to 
extract their generic values. Following [7, 22], we use a discrete algorithm developed by [23] in 
order to compute the sound-channel quasinormal frequencies of the coupled electromagnetic 
and gravitational perturbations of the extremal RN-AdS4 black hole. For numerical purposes 
it is convenient to switch to a new radial coordinate z = 1/u. In terms of this coordinate the 
black hole horizon is located at z = 1, while the asymptotic boundary is located at z = 0. 
The equations (4.20) then read 

z^ [f{z)^'i{z)]' - z^J±{z)-^'^{z) + i^±(z)*±(z) = , (5.1) 

where the prime now denotes derivatives with respect to z, and f{z) is given by 

fiz) = l-{l + Q^)z^ + Qh^. (5.2) 

Similarly, K±{z) and J±{z) can be read off from (4.21)-(4.22). In the method of [23], one first 
factorizes the leading behaviors of the normalizable modes of ^±{z) at the horizon and the 
boundary (which in practice amounts to imposing asymptotic Dirichlet boundary conditions). 
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What remains is then approximated by a power series in z around a point zqi chosen in such 
a way that the radius of convergence of the series reaches both the horizon and the boundary. 
In our case we can take zq = 2- Thus, we write 

£0 r ^ / 1 \ "^ 

^^(^)=e^yIIa^/(,)-^^^a±(^,q) z-;^ . (5.3) 

m=0 ^ ^^ 

The first two prefactors are chosen to select the infahing wave at the extremal black hole 
horizon, as it can be easily seen by rewriting them in terms of infalling Eddington coordinates. 
Substituting (5.3) into (5.1), we obtain a set of M + 1 linear equations for M + 1 unknowns 
{ap(tx>,q)}, that can be written in matrix form as 

M 

^^±p(m,q)aj(nj,q) = 0. (5.4) 

m=0 

Suppose now we choose a specific value for the spatial momentum q. The quasinormal 
frequencies t»^ for that particular value of q are then solutions of 

det ^=^(t»^,q) = 0. (5.5) 

Since the fields ^-t in (5.3) are approximated by a power series, the more terms kept in 
the series, the more accurate the results obtained for the quasinormal frequencies. The 
quasinormal spectra for q = 0.5 obtained as described above is shown in Figure 1 below. 

A distinctive feature of the spectrum at zero temperature is that the series of purely 
imaginary frequencies on the negative imaginary axis, which are evenly spaced at finite tem- 
perature,^ now coalesce to form a branch cut. Indeed, this was observed in the correlators of 
conserved currents in the shear channel [7] as well as in those of chargeless scalar operators 
[22]. As discussed in [7], at zero spatial momentum (q = 0) this feature is explained by the 
explicit appearance of logarithmic terms of the form t»"logtx) (with n G Z) in the Green's 
functions, while for q 7^ it is due to terms of the form t»^'^+ (c.f. (4.72) and (4.73)), where 
2vj^ is generically an irrational number. On the gravity side, one can understand the appear- 
ance of the branch cut by using the tortoise coordinate to write the equations for ^-t in a 
Schrodinger-like form. Then, the non-analyticity is due to the fact that the effective potential 
near the horizon decays as a power law rather than exponentially (as would be the case at 
finite T) , a typical feature of extremal black holes (see [24] and references therein) . That we 
can still observe a discrete series of poles as we move down the imaginary axis is an artifact 
of the numerical algorithm. As we increase the size M of the matrix A defined in (5.4), these 
poles move upwards and become dense over a larger set of imaginary frequencies. The plots 
show only the lower half of the complex-tx) plane, because no instabilities exist. This was 



This has been determined by an explicit calculation of the quasinormal frequencies at finite temperature, 
which for the sake of brevity we do not display here. 
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Figure 1: Sound-channel electromagnetic and gravitational quasinormal frequencies for the extremal 
Reissner-Nordstrom AdS4 black hole. Plot (a) shows the quasinormal frequencies of ^+, with the 
sound modes depicted in blue. Plot (b) shows the quasinormal frequencies of ^_. For both plots, 
q = 0.5 and M = 300. Plot (c) shows the spectra of ^+ and ^_ superimposed. 



also the case for the shear-channel modes, where moreover it was relatively easy to prove the 
stability analytically, at least over a broad range of momenta [7] . 

Besides the coalescent poles giving rise to the branch cut, we observe a discrete series of 
quasinormal frequencies with a non- vanishing real part. Notice that they appear in conjugate 
pairs with the same imaginary part, but with opposite sign for the real part. This was to be 
expected due to the unbroken parity symmetry of the dual theory, a feature that would be 
absent in the presence of a magnetic field, for example. In figure 1 (a) we can distinctively 
appreciate the lowest lying pair (depicted in blue). As we will explicitly confirm below by 
analyzing their dispersion relation, this pair of poles correspond to the sound modes (in the 
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sense described in section 2.2: they satisfy uj{k — )• 0) = ±Csk — iYgk'^ and dominate the decay 
of the longitudinal fluctuations over long wavelengths and time scales). While the imaginary 
part of these quasinormal frequencies is non-zero, it cannot be appreciated in the figures 
above as a consequence of the rather small attenuation constant Tg we will find below. The 
rest of the QNMs with a non- vanishing real part are the so-called overtones, whose dispersion 
relation will seen to be quite different from that of the sound modes. They correspond to 
higher resonances whose effects we expect to become noticeable as we go beyond the small 
frequency and momentum approximation. As it can be observed in figure 1 (c), the overtones 
of ^-|_ and ^_ coalesce as we move down in the complex frequency plane. 

5.2 Dispersion relations: speed of sound and sound attenuation 

At finite temperature, we can easily compute the speed of sound Cg in the boundary field 
theory by using the one-point function of the energy-momentum tensor operator (T^i/). The 
only non- vanishing components of the second fundamental form in our background are 

A-',^4<^^, K%-K',--^, (5.0) 

and with this result we can evaluate (4.34) to obtain 

{T\)^-e = -^^M, (5.7) 

(n> = (n)^^ = i = ^^' (5-8) 

where e is the energy density that had been introduced in (3.6) and P = e/2 is the pres- 
sure. Notice that the expectation value of the trace of the stress-energy tensor {T\) is zero. 
Substituting the values of e and P into the thermodynamic relation c^ = dP/de we obtain 

0^ = ^=^ (59) 

which is the expected result in two spatial dimensions. Note, however, that this result is 
independent of whether the RN-AdS4 black hole background is extremal or non-extremal. 
Thus, for our zero-temperature dual field theory at finite U(l) charge density, the speed of 
sound should also be l/v2. On the other hand, calculating the QNM spectrum for various 
values of q we can numerically construct the dispersion relation of the sound modes. The 
result is shown in figure 2. We can observe that both the real and imaginary parts of vOg 
approach zero as q — )■ 0. Moreover, for small values of the momenta the dispersion relation has 
the form (2.5); fitting the numeric data for q < 0.25, we obtain Cg ~ 0.704 and fiTg ~ 0.083. 
Notice the agreement between the numerically obtained value for the speed of sound and the 
value Cg = l/\/2 ~ 0.707 predicted by (the zero-temperature limit of) hydrodynamics. To 
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Figure 2: Dispersion relation of the sound modes. Plot (a) depicts the imaginary part of rOs as a function of 
momenta, while plot (b) shows the real part of t»s- The dots represent the actual values obtained numerically, 
while the solid line results from fitting the data by a smooth curve. The inset plots compare the expected 
dispersion relation (5.12) (red line) with the values obtained numerically (black line), for small values of the 
momenta. For large momentum, the slope of the real part becomes unity within numerical precision. 



gain some insight about the value obtained for the sound attenuation constant Tg, we can 

again investigate the T — )• hniit of weh- known hydrodynamical relations. As discussed in 

section 2.2, such a procedure may or may not yield the right value for the theory at zero 

temperature. Putting aside this issue for the moment, using the ratio s/e = 2ttQ / {^M) (see 

(3.6)) we find 

/i r] Irj fis 2Tr T] Q 

2e + P ~3sT~YsM" 



AtT^ 



(5.10) 



At extremality, Q^ = 3, M = 4, and the ratio of the shear viscosity rj to entropy density s 
was found to be r]/s = l/(47r) [6], just as in the well-known finite temperature examples (see 
[8] and references therein). Thus, by taking the T — >• limit of (5.10) we expect 



fiV 



1 



«lT->0 



8\/3 



0.072. 



(5.11) 



Summarizing, by taking the zero temperature limit of hydrodynamics we anticipate the sound 
modes (the one with positive real part, say) to have a dispersion relation of the form 



t»^(q -> 0) 



V2 



8^3 



(5.12) 



The inset plots in figure 2 compare the dispersion relation (5.12) (red line) with the one 
obtained from the actual numerical data at T = (black line). While the agreement in the 
real part of the frequency prevails in a relatively broad range of momenta, the imaginary 
part computed by the different methods agrees only in the low frequency and momentum 
regime. In particular, the actual value fiTg ~ 0.083 for the sound attenuation constant differs 
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slightly from the expected value (5.11). This discrepancy might be attributed to numerical 
uncertainties, but we note that it is rather large, of order 13%, in fact much larger than the 
uncertainty implied by the results for the speed of sound. 

Turning now to the overtones, we can also follow their trajectories in the complex fre- 
quency plane as we vary the momentum q and determine in this way their dispersion relation. 
Figure 3 shows the result for the first five overtones in the spectrum of ^+. 
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Figure 3: Dispersion relation of the first five overtones of \['+. Plot (a) shows the variation of the imag- 
inary part of the corresponding quasinormal frequencies as we vary the momentum q. Plot (b) shows the 
corresponding change in the real part of the frequency. 

As mentioned above, these frequencies do not approach zero as q — )• 0, so these modes are 
not "hydrodynamic" in the broad sense. In particular, notice that the imaginary part of 
the quasinormal frequencies is approximately constant for small values of the momenta with 
q < 1 (fc < /u), and also for large values of the momenta of the order q ~ 10. In contrast, 
the real part of these overtones is a monotonically increasing function of q , satisfying a linear 
dispersion with unit slope (within the numerical precision) for values of the momenta greater 
than q ~ 7. This feature is explained by the fact that at large q and m the influence 
of the charge density becomes less important, and the background effectively approaches 
a Lorentz-invariant vacuum. Provided that the higher resonances behave like sharp, well- 
defined particle-like excitations at large momenta, their dispersion relation becomes effectively 
relativistic in the UV. The same qualitative features are observed for the first five overtones 
of ^_, whose dispersion relations are shown in figure 4. 

5.3 Residues 

In order to further understand the field theory spectrum, it is important to compute the 
residues of the retarded correlators at the various quasinormal frequencies, as they determine 
the weight with which they contribute to the spectral function. To compute the residues 
numerically, we found it convenient to use a method combining series expansions and numeric 
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Figure 4: Dispersion relation of the first five overtones of ^_. Plot (a) shows the variation of the imag- 
inary part of the corresponding quasinormal frequencies as we vary the momentum q. Plot (b) shows the 
corresponding change in the real part of the frequency. 



integration. Basically, the idea is to approximate by a power series the solution for the fields 
^±{z) near the horizon (i.e. for z > 1 — e) and near the boundary (i.e. for z < 6), and 
integrate numerically from z = 6toz = l — e. In this way, one avoids the stiffness issues 
associated with numeric integration near the singularities.^ In order to ensure convergence 
of the series expansions, for the values of the momenta we have studied it was sufficient to 
choose e = 10~^ and 6 = 10~^, and we kept ten terms in the series approximations at each 
end of the integration interval. By matching the numerically obtained solution for the fields 
and their derivatives with the series expansions at z = 6, we can construct the coefficients 



,(0) 



,(1) 



^^-^ = ^-1- and ^^ in (4.23) as a function of the frequency (for fixed q) and with them we 
obtain the canonical momenta II-i- using (4.27). At this step we can compute the residue of 
II-i- at a given quasinormal frequency t»^, which is given by 



Res n-i-(ix)^ 



1 ^ 



(1), 



m 



C±5t»^±(t») 



(5.13) 



In practice, the denominator in this expression can be evaluated by discretizing the to- 
derivative. Prom the residue of II-i- we can calculate Res Gi and Res G2 using (4.39)-(4.40), 
which in turn allows us to compute the residue of the actual retarded Green's functions 
constructed from Gi and G2 as given in (4.41)-(4.53). 

In principle one should not only compute the residues of the sound modes, but those of 
the higher resonances as well. For example, knowledge of the residue at the first overtone 
allows one to compute the relaxation time scale into the hydrodynamic regime [26, 9]. While 

^This method was used in [25] to compute the quasinormal spectrum. In particular, the quasinormal fre- 
quencies are obtained as the zeros of the Wronskian between the numeric solution and the series approximation 
at the matching point. 
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being cognizant of the importance of these residues, due to the complexity of the numerical 
calculations involved, here we will restrict ourselves to the hydrodynamic (sound) modes. As 
a representative of our results, in figure 5 we display the absolute value of the sound pole 
residue for Gu tt and Gx x- 



(a) 



{b) 
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Figure 5: Absolute value of the sound pole residue as a function of momentum. Plot (a) shows the residue 
of Gtt,tt ; plot (b) shows the residue of Gx,x- 

In particular, we observe that the absolute value of the residue approaches zero as q — )■ 0, 
and behaves as a monotonically increasing function of q . 

5.4 Green's functions 

By using the same combination of series approximations and numeric integration that was 
employed in the calculation of the residues, one can numerically construct the retarded Green's 
functions (4.41)-(4.53) as a function of frequency (for fixed q). Since we have also computed 
the sound pole residue, one can compare the Green's functions obtained numerically with 
an approximation of the form (2.6). The results for the imaginary and real parts of Gtt,tt 
at q = 0.5 and q = 1 are shown in figure 6. It is worth mentioning that the contact terms 
presented in appendix B should be taken into account when plotting the real part of the 
retarded correlators. 

At q = 0.5 (plots (a) and (b)) we observe an excellent agreement between the Green's 
function obtained numerically and the approximation obtained by truncating the spectrum 
to the sound modes only. Moreover, the sound pole resonance is sharply peaked, so that 
these modes are well-defined particle-like excitations. This is natural, since at small values 
of q the modes with oj{k — )■ 0) — )■ are expected to dominate the spectrum. For q = 1, 
the real part of the Green's function is also accurately described by the sound mode. In 
contrast, the numerically computed imaginary part shows two smaller peaks between tx) = 1.0 
and vo = 1.2, which are not captured by the truncated approximation including the sound 
mode only (dashed line). As we increase the momentum, we should in general expect the 
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Figure 6: Real and imaginary parts of the Gtt, tt Green's functions as a function of the (real part of the) 
frequency. Plots (a) and (b) correspond to q = 0.5, while plots (c) and (d) correspond to C) = 1. The dots 
represent the actual numerical data, and the solid (red) is the interpolation of the data by a smooth curve. 
The dashed (black) line represents the Green's function approximated by (2.6), with the spectrum truncated 
to include the sound modes only. 

contribution of the higher resonances to become noticeable, and it is therefore tempting to 
attribute these peaks to the effects of the first overtones of ^4. and ^_ , respectively. However, 
the position of the peaks in the plot slightly differs from the first overtone frequencies as 
computed using Leaver's method. This might indicate that the numeric precision in our 
computations might not be adequate to conclusively identify the effects of higher resonances. 
A possible way of elucidate the nature of these peaks is to include the contribution of the first 
overtones in the spectral approximation (2.6) and compare it with the curve obtained from 
the numerics. Naturally, this approach presumes a knowledge of the residues of said poles 
with high enough accuracy, which we were not able to obtain with the numeric techniques 
described above. 
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6. Conclusions 

Using gauge/gravity duality techniques, we have constructed the correlation functions of the 
conserved charge current and energy-momentum tensor of a (2-|-l)-dimensional strongly cou- 
pled field theory at finite density and zero temperature. This problem is more challenging 
than the cases addressed previously in the literature: not only we have considered the coupled 
electromagnetic and gravitational fluctuations, but we also faced the additional complications 
introduced by the more severe singularity structure of the equations of motion in the extremal 
background. In the process of these numerical computations, we resolved a number of techni- 
cal obstacles, and in particular were able to set up the Dirichlet problem for decoupled master 
fields that determine the correlation functions of interest. We have found that the correlation 
functions at finite density and zero temperature are essentially hydrodynamic in character 
at low frequency and momentum, with the chemical potential setting the scale (rather than 
temperature). In the IR, they satisfy scaling behavior originating from the IR CFT dual to 
the near-horizon geometry, and at generic momenta display a dominant sound mode and a 
sequence of resonances (with typically small residue), along with a branch cut at tj = 0. It 
would be interesting to apply these techniques to other systems of direct physical interest. 
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A. Master field equations 

In this section we provide further details about the master fields introduced in the body of 
the paper. 

A.l Charged black hole 

The Einstein-Maxwell equations (3.15)-(3.24) can be reduced to two decoupled second-order 
differential equations for the so-called master fields ^±{u) [11]. In order to obtain the master 
fields, define first 

A{u) = v?a[{u) + iihyy{u) fih\{u) , (A.l) 

^{u) = uhyyiu) - H{u)-\^f{u) [/i^;(n) + /i%(n)] , (A.2) 
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where H{u) = v? f'{u) + Q^<\^. Equations (3.15)-(3.24) now reduce to a set of two coupled 
second-order differential equations 

[u^f{u)A'iu)]' + 2fif{u)<^'{u) + Ci{u)A{u) + fiC2iu)<^{u) = 0, (A.3) 

[u^f{u)^'{u)]' + C3(n)$(n) + fi-^Q^C4iu)A{u) = 0, (A.4) 



with 








(A.5) 


u^ J 


(A.6) 




(A.7) 


A-^t). 


(A.8) 



v?f{u) \ ■" ' H{u) 

C2{u) = -^ {[QW - 10n2/(^)] f[u) + 24n/(n)[l - /(n)] + 

Ca{u) = -j^ (^n^fiu) [Wf{u) + nf\u)] + 24n/(u) [f{u) - 1] 

where 

V{u) = -^±-^[9u''f'{uf + u' [3(22/(n) - 8^ + Q^ (9q2^2 ^ 4)] ^/(^)2 

+ n3 [216n6/(n)(/(^) - 1) + 8Q2c,2^4(2/(^) _ 3) + ^4^2 (^2^2 ^ g^] ^/^^^ 

+ 288nV(n)(/(n) - l)^ - 2AQWu^f{u){f{u) - 1) + Q'^q^ (q^n^ + 4) }. (A.9) 

Equations (A.3) and (A.4) can then be decoupled by introducing the master fields 

$±(-u) = Q±(u)$(u) + fi-^Q^A{u) , (A.IO) 

where 

a±{u) = ^F±(Q) - ^, F±(g) = ^ [(1 + Q2) ± 7(l + Q2)2 + (i6/9)Q4q2j . (a.h) 

To obtain the equations for the master fields, solve (A.IO) for A{u) and <I>(u) in terms of 
^± (u) to obtain 

H-'^Q'^A{u) = 7 [a+(u)$„(n) - a_(u)$+(u)] , (A.12) 

$(-u) = 7 [$+(n) - $_ (u)] , (A.13) 

where 

1 /I 

(A.14) 



a+{u)-a-{u) 3v'(l + Q2)2 + (16/9)Q4q2 
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Substituting (A. 12) and (A. 13) into (A. 3) and (A. 4), we arrive at 

[u^f{u)'^'±{u)]' ± 7 [Q^f'iu) + a±{u)Ci{u) - Q^C2{u) - a^Csiu) + al{u)C^{u)] c^^(n) 

T 7 [Q'^f'iu) + a^{u)Ci{u) - Q^C2iu) - a±C3{u) + a+(n)a_(n)C4(tx)] $±(u) = . 

(A.15) 

Note that 

Q^f'{u) + a±{u)Ci{u) - Q^C2{u) - a±{u)Cz{u) + a\{u)Ci{u) = . (A.16) 

Simplifying the equations for the master fields, we arrive at 

u^f{u) [u'f{u)^'^{u)\ + [Q2«,2 _ u^{u)] $±(n) = 0, (A.17) 

where 

U±{u) = ±-ff{u) {a±{u)V{u) + -\- [QWH{u) + 8QV(n)] + — [H{u) - 2QW] 

^H^) [48«'/(^)(/H - 1) + 20txV(n)/'(n)] | . (A.18) 

A. 2 Chargeless case 

In the case where the background in chargeless, i.e. Q = 0, the gravitational and electro- 
magnetic perturbations decouple, and the fields $± reduce to a single master field. In our 
notation, this field is given by <I*+(u) = <I>(n) and it encodes the information from the gravity 
sector only. It can be easily seen to satisfy 

u^f{u) [u^f{u)^'{u)^ + [9-^2 - U{u)\ $(n) = , (A.19) 



where 



f{u) 



U{u) = ^-^ 2 [2«^/H + 3 (1 + 3g2u) (1 + 9qW)] , (A.20) 



and 



w=——, q = -r^- A.21 



B. Contact terms 



Here we display explicitly the contact terms that supplement the Green's functions (4.41)- 
(4.53). For notational simplicity, in the expressions below we suppress the overall constant C 
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coming from the normalization of the boundary action. We display the results for arbitrary 
temperature (equivalently, Q), but it is understood that in the calculations of the paper we 
have used the extremal value Q = \/3. The contact term contributions to the {T^uTpa) 
correlators are-*^*^ 



(1 + Q^) 2tp4 - 7q^to2 + 8q^ _ (1 + Q^) m* - 2q'^m'^ - 2q^ 

2 (t»2-q2)2 ' ^^^,tt- ^ (u)2-q2)2 

(1 + Q^) P)^ + 2q=^ 3(l + Q^)qto(to^-2q=^) 

2 t»2_q2 ' <~^xt,tt- 2 (a)2-q2)2 

(l + Q^)p)^-5q=^tX)^+q^ _ (]^^2^^+q2 

2 (tx>2-q2)2 ' ^y?;,xx- 2 «>2 _ ^2 ' 



3(1 + Q2) q3tx> _ (1 + Q ^ 

(0,2 - q2)2 ' ^yy^yy - 2 



_ oyx -r ^^ ; vi w _ V^ -r S^ ; /p .X 

'-'xtjxa; — r, /.„5 2A2 ' ^yy,yy ~ o ' vJ-"-^^ 



3(1 + Q^) q»x> ^ _n^n2^2^!-4^'«''-^' 



<--a;t,j/y - 7) ™,2 _ „2 ' '-'xt.xt - 1,1 + V j 7Z;2Tr;;2\2 ' ^^•^'' 



On the other hand, there is no contact term contribution to the current-current correlators 
of the form (J^Jjy): 

Ci,t(tx>,q) = 0, C^,i(t»,q) = 0, C^,^(t»,q) = . (B.6) 

For the "mixed" correlators {JpTup) we find 

_ 4Q2q>^^q2) _ 4Q2 q^(^^^q2) 

'** ~ ^ (-Q,2_q2)2 ' '^x,tt-—. 7172 l2\2~' V^-'^ 

_ 4Q^ q^ __^__^ 

"""i'^^ ~ „ (jp2_q2)2' L-x.xa; - ,. ^_2 ^2^2 ' V^' 



_4Q2 qn,(»,2-2q2) _ .^ „ ,„ ■. , 

^*'"*"^^ (tx,2-q2)2 ' C,,xt--— (j„2_q2)2 • i^-lU) 



4Q2 


(t»2 -q- 
q^tx) 


2)2 ' 


4Q2 


(t»2 -q: 
qt» 


2\2 ' 


;U t»2 — q2 ' 
4Q2tx,2(n,2_ 


-2q2) 



^°The name "contact term" usually denotes contributions which are analytic in frequency and momentum. 
As noticed in the body of the paper, the apparent pole at t» = zbq is spurious and denotes a choice of 
normalization. In this sense, the expressions presented in this appendix are to be regarded as true (analytic) 
contact terms. 
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